The effects of Mn concentration on spin-polarized transport 
through ZnSe/ZnMnSe/ZnSe heterostructures 

Alireza SaffarzadehQ 

Department of Physics, Tehran Payame Noor University, 
Fallahpour St., Nejatollahi St., Tehran, Iran 
(Dated: February 6, 2008) 

Abstract 

We have studied the effects of Mn concentration on the ballistic spin-polarized transport through 
diluted magnetic semiconductor heterostructures with a single paramagnetic layer. Using a fitted 
function for zero-field conduction band offset based on the experimental data, we found that the 
spin current densities strongly depend on the Mn concentration. The magnitude as well as the 
sign of the electron-spin polarization and the tunnel magnetoresistance can be tuned by varying 
the Mn concentration, the width of the paramagnetic layer, and the external magnetic field. By an 
appropriate choice of the Mn concentration and the width of the paramagnetic layer, the degree of 
spin polarization for the output current can reach 100% and the device can be used as a spin filter. 
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I. INTRODUCTION 



Recent advances in controlling the spin-polarized transport in metallic magnetic struc- 
tures, and successful applications of such systems in electronic devices renewed world- 
wide interest in diluted magnetic semiconductor (DMS) structures. These materials are 
a promising component for the new spin-based information technology, in which the spin 
degree of freedom of the electron can be utilized to sense, store, process and transfer infor- 
mation. II- VI DMSs 0,0], are a class of materials that has received much attention in recent 
years due to their interesting physical properties and the potential applications in integrated 
magneto-optoelectronic devices P] . They are known to be good candidates for effective spin 
injection into a nonmagnetic semiconductor (NMS) because their spin polarization is nearly 
100% and their conductivity is comparable to that of typical NMS. Among different types of 
II- VI DMSs, (Zn,Mn)Se is a very promising one for spin injection, which has been previously 
used for spin injection experiments into GaAs [Q], and ZnSe 

HQS. 

All A^L^jMuajB^^ alloys are direct-gap semiconductors like their A^^B^^ parent materials. 
In the presence of a magnetic field, the band edges in these materials undergo a huge spin 
splitting due to the sp — d exchange between carriers and localized magnetic ions ^, while 
the splitting in the nonmagnetic layers is much smaller. When the large Zeeman splitting in 
the magnetic layers overcomes the band offsets in both conduction and valance bands, the 
magnetic layers act as barriers for electrons and holes in the spin-up state, and as quantum 
wells for the spin-down state. 

In recent years, the spin-polarized transport in II- VI DMS heterostructures has been 
investigated theoretically. Based on a quantum theory and the free-electron approxima- 
tion, Egues jlj, Guo et al. and Chang and Peeters 1^ studied spin filtering in 
ZnSe/Zni_j.Mna,Se/ZnSe heterostructures in the ballistic region. The results showed a strong 
suppression for one of the spin components of the current density with increasing the exter- 
nal magnetic field. Also, Zhai et al. [isl investigated the effects of conduction band offset 
on the spin transport in such heterostructures and showed that the positive zero-field band 
offset can drastically increase the spin polarization. However, they have not considered the 
effects of Mn concentration and the width of the paramagnetic layer on the spin-polarized 
transport. Thus other aspects of these heterostructures remain to be explained. 

In the present work, we study theoretically the dependence of spin-polarized transport 
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on Mn concentration in ZnSe/ZnMnSe/ZnSe structures. The ZnSe layers are considered as 
emitter and collector attached to external leads. We assume that the carrier wave vector 
parallel to the interfaces and the carrier spin are conserved in the tunneling process through 
the whole system. These assumptions can be well justified for interfaces between materials 
whose lattice constants are nearly equal; and when the sample dimensions are much smaller 
than the spin coherence length. Using a quantum theory, we study the effects of Mn con- 
centration and the width of the paramagnetic layer on spin-dependent current densities, the 
degree of spin polarization, and the magnetoresistance ratio. 

In section 2, we present a fitted function for zero-field conduction band offset and describe 
the model. Then, the spin current densities, the electron-spin polarization, and the magne- 
toresistance ratio for ZnSe/ZnMnSe/ZnSc hctcrostructurcs arc formulated. In section 3, the 
numerical results for the above quantities are discussed in terms of the Mn concentration, 
the width of the paramagnetic layer, and the applied voltage. The results of this work are 
summarized in section 4. 

II. MODEL AND FORMALISM 

Consider a spin unpolarized electron current injected into ZnSe/Zni_^Mn2.Se/ZnSe struc- 
tures in the presence of magnetic and electric fields along the growth direction (taken as z 
axis). In Mn-based DMS systems, the conduction electrons that contribute to the electric 
currents, interact with the 3d^ electrons of the Mn ions with spin -S" = | via the sp-d exchange 
interaction. Due to the sp-d exchange interaction, the external magnetic field gives rise to 
the spin splitting of the conduction band states in the Zni_3;Mna;Se layer. Therefore, the 
injected electrons see a spin-dependent potential. In the framework of the parabolic-band 
effective mass approximation, the one-electron Hamiltonian of such system can be written 
as 

1 eV 7 

^ ^ + + + ^^^^^ + ^-^^^ " ^ ' 

where the electron effective mass m* is assumed to be identical in all the layers, and the 
vector potential is taken as A = (0,i?a;, 0). Here, Vg = ^gsl^-BC ■ B describes the Zeeman 
splitting of the conduction electrons, where a is the conventional Pauli spin operator; Vx{z) 
is the heterostructure potential or the conduction band offset in the absence of a magnetic 
field, which depends on the Mn concentration x and is the difference between the conduction 
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band edge of the Zni_a;Mna;Se layer and that of the ZnSe layer; Va^{z) is the sp-d exchange 
interaction between the spin of injected electron and the spin of Mn ions and can be cal- 
culated within the mean field approximation. Hence, the sum of the last two terms can be 
written as 

V.{z) + V^^{z) = [{l/2)AE{x) - Noaa,Xeff{S,)]e{z)e{d - z) , (2) 

where 

AE(x) = -0.63x + 22x^ - 195x^ + 645x^ . (3) 

Here, AE{x){= Eg{x) — Eg{0)) is the sum of the conduction and valance band offset under 
zero magnetic field, when the real Mn concentration is x. In Fig. 1, we have shown the Mn 
concentration dependence of the energy gap Eg{x) in Zni_^.Mn^Se layer. It is clear that the 
band gap of Zni_a;Mn2;Se varies anomalously with x; it shows a minimum at x ~ 0.02 and 
increases linearly with x for x > 0.05. This anomalous behavior of Zni-^jMn^^Se which called 
band-gap bowing, probably refers to the sp-d exchange interaction y. In deriving AE{x), 
we have used of the fitted curve of Eg{x). We have approximated AE{x)/2 as the zero-field 
conduction band offset. This function does not depend on spin, and for x < 0.043 behaves 
as a quantum well, and as a potential barrier for x > 0.043. 

In Eq. (121), A^oa is the electronic sp-d exchange constant, = ±1/2 (or t, |) are the 
electron spin components along the magnetic field, Xeff = x{l — x)^^ is the effective Mn 
concentration used to account for antiferromagnetic Mn-Mn coupling, and Q{z) is the step 
function. (Sz) is the thermal average of zth component of Mn^+ spin in the paramagnetic 
layer which is given by the modified | Brillouin function SBs^^bB /ksiT + Tq)], where Tq 
is an empirical parameter representing antiferromagnetic interactions between the Mn ions 
jlS^ . The last term in Eq. (|H) denotes the effect of an applied voltage Va along the z axis 
on the system, and d is the width of the paramagnetic layer. It is important to note that, d 
is much smaller than the spin coherence length in the semiconductors. Therefore, we have 
neglected the effects of spin-flip processes in the Hamiltonian of the system. 

In the absence of any kind of scattering center for the electrons, the motion along the 
z-axis is decoupled from that of the x-y plane. Therefore, in the presence of magnetic field 
S, the in-plane motion is quantized in Landau levels with energies En = {n + l/2)hujc, where 
n = 0, 1, 2, ■ ■ ■ and Uc = eB/m*. In such case, the motion of electrons along the z axis can 
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be reduced to the following one- dimensional (ID) Schrodinger equation 

- 2^%^ + = W , (4) 

where is the longitudinal energy of electrons and Uo-^ {z, B) = Vs + Vx{z) + V„.^ {z) — eVaz/d 
is the effective potential seen by a traverse electron, which includes the effects of spin, 
conduction band offset, and external magnetic and electric fields. The general solution to 
the above Schrodinger equation is as follows: 



(5) 



where = ^2m*{E, - K)/^, ha, = ^'2m*{E, - K + eK)/^; Ai[p,^(2;)] and Bi[p^^(2;)] 
are Airy functions with p^,(z) = [d/ {eVaX)]{E^ - U^J, and A = [-h^d/ {2m*eVa)Y^^] ^j^z 
and Bj(j^ (with j=l-3) are constants which can be obtained by matching the wave functions 
and their derivatives at the interfaces of Zni-xMn^^Se and ZnSe layers. The matching results 
in a system of equations, which can be represented in a matrix form jlfl] . 




1 (0) M2 (0) M2~i (rf) M3 (rf) 



A 



3cTz 



total 




(6) 



Here, Mtotai is the transfer matrix that connects the incidence and transmission amplitudes, 
and 



M,{zi) 



i)Uz) (z) 



(7) 



\ dz 



dz 



where, ipfiz) and ipj{z) are, respectively, the first and second term of the wave functions 
in each layer, without considering their coefficients. Therefore, the transmission coefficient 
of the spin electron, which is defined as the ratio of the transmitted flux in the collector 
to the incident flux in the emitter, can be written as 

2 



T^^{E,,B,Va) 



hoz 



1 



^^■^ total 



(8) 
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where M^j^^i is the (1,1) element of the matrix Mtotai- The spin-dependent current density 
can be determined by 



°o r+oo 

J^XB) = JoBj2 T^AE.,B,Va) 



x{f[E, + (n + ^)hu, + Vs] - f[E, + {n + ^)hu, + Vs + eVa]}dE, , (9) 

where Jq = e^/47r^^^, f{E) = 1/{1 + exp[(£' — Ep) /ksT]} is the Fermi-Dirac distribution 
function in which is the Boltzmann constant, T is the temperature, and Ep denotes the 
emitter Fermi energy. 

The degree of spin polarization for electrons traversing the heterostructure is defined by 

^ MB) -MB) 

MB) + MB) • ^ ' 

where J-^ ( J|) is the spin-up (spin-down) current density. On the other hand, in the absence 
of external magnetic field, the Zeeman splitting of the conduction electrons Vg and the 
spin dependent potential Va^{z) are zero. In this case, the effective potential reduces to 
Ua^{z,0) = Vx{z) — eVaz/d which is spin-independent. Thus, by using Ua^{z,0) and a 
procedure completely analogous to the one used for the case of i? 7^ 0, one can obtain the 
following formula for the total electric current density 



ATT^h^ J Jo " ' \l + exp[{EF- E,-eVa)/kBT] 

The factor two is due to the degeneracy of the electron spin in the case of 5 = 0. Although 
the transverse momentum ky was not appeared in the above equation, the effects of the 
summation over ky have been considered in our calculations. Here, we mention again that 
the effective mass is independent of layer. When taking both the transverse motion and 
the layer-dependent effective mass of the electron into account, the transmission coefficient 
can have a pronounced dispersion in ky space In this case, one cannot simply reduce 
the 3D Schrodinger equation to the ID equation and integrate the ky to obtain the current 
density, as we did here. 

The linear conductances per unit area are given by G{0) = J{0)/Va and G{B) = 
J2a^ JaAB)/Va for 5 = and B ^ 0, respectively. Therefore, the tunnel magnetoresistance 
(TMR) or magnetoresistance ratio in such heterostructures can be described quantitatively 
by 

G(0)-G(B) 
= G(B) 



= (12) 

In next section, we will present the numerical results for J^., P and TMR in terms of Mn 
concentration. 



III. RESULTS AND DISCUSSION 

In our numerical calculation we have taken the following values: m* = 0.16 me {nif, is 
the mass of the free electron), Ep = 5 meV, B = 4 T, Qs = 1.1, T = 2.2 K, To = 1.7 K, and 
Noa = —0.26 eV [l^. Fig. 2(a)-(d) show the dependence of spin-up current densities on the 
Mn concentration for several choices of d and V^. In all of the widths of the paramagnetic 
layer, for a fixed value of Mn concentration, the spin-up current density increases with the 
applied voltage. With increasing the Mn concentration [x > 0.043), the spin-up electrons 
see a higher potential barrier, while with increasing the width of the paramagnetic layer, 
these electrons see a thicker potential barrier. Thus in both cases, the tunneling probability 
for this group of electrons decreases and this leads to a reduction of the current density. 
From the figures, we also find that, for d =200 nm and x > 0.055, the transmission for 
spin-up electrons is completely suppressed. On the other hand, with increasing the applied 
voltage, the potential barrier tilts and the effective width of the barrier becomes narrower; 
therefore, the spin-up current density increases, as shown in Fig. 2. 

The dependence of spin-down current densities on the Mn concentration is shown in Fig. 
3(a)-(d) for several values of d and Va. For small widths of the paramagnetic layer, the 
variations of the current density are relatively small with an oscillatory behavior. These 
variations increase with the width d. Qualitatively, the voltage dependence of spin-down 
current density is similar to the spin-up one. This means that both spin current densities 
increase with the applied voltage. We should note that at x = the current densities are 
nearly sp in- independent , because is zero and Vg is very small. With increasing d and x, 
peaks are observed in the spin-down current densities. The reason is that, the paramagnetic 
layer behaves as a quantum well for spin-down electrons; thus, the enhancement of the width 
of the paramagnetic layer, varies and shifts the position of the resonant states formed in 
the well to the lower energy region. This leads to the formation of peaks in the spin-down 
current densities l20|. 
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In Fig. 4(a)-(d) we show the spin polarization for electrons traversing the heterostructure 
as a function of x for various values d. According to the above discussion, for x = the 
spin-up and spin-down current densities have equal values, hence the spin polarization is 
zero. As the Mn concentration increases, the spin-down current density first decreases, while 
the spin-up one is approximately constant; therefore, the spin polarization is negative. On 
further increasing the Mn concentration, the spin-up current density decreases exponentially; 
however the spin-down one can increase and when J| vanishes, has non-zero values. 
Consequently, the spin polarization becomes positive and increases for these values of the 
Mn concentration. It is also important to note that, for widths near 200 nm, the spin 
polarization approaches 100%, which shows that, the system acts as a spin filter. 

The magnetoresistance ratio of ZnSe/Zni_a,Mna.Se/ZnSe heterostructures is another 
quantity which is sensitive to the Mn concentration. In Fig. 5(a)-(d) the TMR is plot- 
ted as a function of x. With increasing the width of Zni_2.Mna,Se layer, a peak appears in 
the TMR curves. This peak is shifted slowly towards lower values of the Mn concentration 
X, and its height increases with the width. From the figures, it is clear that the sign of TMR 
can be positive when G{0) > G{B) or negative when G(0) < G{B). The reason of negative 
TMR can be understood by considering the effects of Mn concentration on spin transport 
in such heterostructures. With increasing the Mn content (x > 0.043), the conduction band 
offset increases and the paramagnetic layer acts as a barrier for both spin-up and spin-down 
electrons in zero magnetic field. On applying a magnetic field, however, the spin-down elec- 
trons see a quantum well (when Vx + V^^ < 0) and the current density for such electrons 
increases. Hence, G'(O) < G{B) and TMR becomes negative. The results also show that, for 
a fixed width of the paramagnetic layer, the electron-spin polarization and the TMR curves 
do not depend strongly on the applied voltage, as shown in Figs. 4 and 5. 

Guo et al. studied the effects of zero- field conduction band offset on spin transport 
with Vx = —5, 0, +5 meV and x = 0.05. They found that, for = +5 meV the spin currents 
are highly polarized, while, for Vx = —5 meV the spin polarization is very low. Our present 
results based on the experimental data, however, indicate that the conduction band offset 
depends on the Mn concentration. Thus, in order to understand the correct description of 
the effects of the conduction band offset on spin transport, the dependence of this quantity 
on the Mn concentration was included. It is clear that the width of the paramagnetic layer 
is also one of the main factors in spin-polarization of the output current of the system. 
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Therefore, the obtained results clearly illustrate that the current densities and hence the 
degree of spin polarization and the TMR can be tuned by changing the Mn concentration 
and/or the width of the paramagnetic layer. 

IV. SUMMARY 

In this paper, using the transfer matrix method and the effective-mass approximation we 
investigated the ballistic spin-polarized transport in ZnSe/Zni_j;Mn-rSe/ZnSe heterostruc- 
tures. We examined the effects of Mn concentration, the width of the paramagnetic layer, 
and the external magnetic and electric fields on the spin current densities, the electron-spin 
polarization and the TMR. The numerical results show that the zero-field conduction band 
offset which varies with the Mn concentration, plays an important role in the spin current 
densities. The spin polarization and the TMR are not very sensitive to the applied voltage. 
However, by adjusting the Mn concentration and the width of the paramagnetic layer, the 
output current exhibits a nearly 100% spin polarization, and also the sign of the TMR can be 
positive or negative. The presented results may be helpful from a technological application 
point of view such as the generation of spin-polarized injection electrons. 
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FIG. 1: Energy gap of Zni-^Miii^Se as a function of Mn concentration x at T=2.2 K. The experi- 
mental data (full squares) is taken from The solid curve is a fit to the data. 
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FIG. 2: Spin-up current densities as a function of Mn concentration x for different applied voltages 
and widths of Zni_j;Mn-rSe layer. 
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FIG. 3: Spin-down current densities as a function of Mn concentration x for different applied 
voltages and widths of Zni_-rMnj;Se layer. 
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FIG. 4: Electron-spin polarization P as a function of Mn concentration x for different applied 
voltages and widths of Zni_3;Mn-rSe layer. 
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FIG. 5: TMR as a function of Mn concentration x for different applied voltages and widths of 
Zni_-rMn3;Se layer. 
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